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(N : An effective action is obtained of a Bose gas in the bulk separated into two regions by a strong 

^ ' external potential depending on the single coordinate. The main attention is focused on the 

^ ' relaxation of the difference between phases of the weakly coupling condensates of the different 

Qbulk domains separated from each other by the external potential. The cases of low and high 
, temperatures are considered. 

in 

(N 

The experimental realization of the Bose-Einstein condensation at ultralow temperatures in atomic vapors [1] 
provides an example of the systems in which the approximation of the weakly non-ideal gas is well applicable. This 
O fact is connected with the small density of particles in the systems concerned. The realization of such systems gives the 
^ possibility of the experimental investigations of the macroscopic quantum phcnomcnons of different types. Recently, 
Ch ' the manifestations of the macroscopic quantum phase and its behavior are studied actively [2-5], [8]. 

The investigation of the kinetic phenomena due to the relaxation of the order parameter is of a doubtless interest 
for the study of macroscopic quantum phenomena. In the present work we study the spatial relaxation of the phase 



^ of the order parameter in the inhomogeneous Bose gas with the weak coupling between different spatial regions due 



to a barrier tunneling. 

The system concerned in the present work is a Bose gas of the small density in the bulk separated into two regions 
by a strong external potential U {'r') which depends on the x coordinate and does not depend on T*^ — {y,z)^ 
'"^ , i.e., U (T*) = U {x). For simplicity, we assume that the external potential has a rectangular shape U {x) — Uq for 
^ ■ —d < X < d and U (x) — beyond this region of the x coordinate. The height of the external potential Uq is supposed 
O J to be the largest energy parameter in the system, in particular, Uq >> n where /i is the chemical potential /i — nX, 
A is the scattering amplitude of Bose particles, and n is the density of the Bose gas. Due to this assumption the 
interaction between particles of the Bose gas in the region of the influence of the external potential, i.e., in the region 
—d < X < d, can be neglected. 

The left and the right domains of the bulk are supposed to have the same temperature T smaller than the Bose- 
■ condensation temperature Tc. We assume that the densities of the right and the left domains of the bulk have the 
values ni = n — ^An and n2 — n + ^Ati correspondingly and the difference between densities An is much smaller 
than the average density n. We will consider both the case of the nonzero An and the case of An ~ 0. At the same 
time, in these two cases the phases of the left and the right Bose condensates are supposed to be different at the initial 
time moment. This assumption means that at the initial time moment there is a nonzero current from one side of the 
bulk to another. This non-equilibrium initial state will relax to the equilibrium state having the same densities and 
the same phases of the condensates for both sides of the bulk if An — or if An can change in time, and will relax 
^ to the stationary state if the density difference An keeps constant and nonzero. This relaxation process is studied in 
I the present work. 

■ O Note that in the case of superconductors the similar initial non-equilibrium state results in the Josephson oscillations 

with a small damping [6] , [7] . The essential distinction of the Bose gas consisting of neutral atoms from superconductors 
in which the Cooper pairs represent the charged objects is the presence of the gapless excitation spectrum. The 
presence of low energy excitations results in the change of the character of the initial state relaxation making the 
relaxation essentially faster. In the case of the homogeneous in y-z plane difference between the phases of the 
condensates the excitations radiated during the relaxation process have the one-dimensional character. The one- 
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^ [ dimensional character of the radiated excitations results in a divergency of the small momentum correlators of these 
- - 1 excitations. Hence, the consideration of the relaxation process to the second order in the tunneling amplitude is not 
correct in contrast to the case of superconductors [7] and requires more accurate analysis [9], [14], [15]. 

I. THE EFFECTIVE ACTION FOR THE DIFFERENCE OF PHASES. 

Below the approach to the problem of the phase difference relaxation without using the perturbation theory in 
the tunneling amplitude is developed. For this purpose the effective action for the difference between phases of the 
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Bosc field at the right-hand and the left-hand boundaries of the potential barrier is derived. This effective action 
is obtained by the integration over the bulk components of the Bose field with the fixed values of this field at the 
potential barrier boundaries. The obtained effective action contains the relaxation part which is proportional to the 
first power of the frequency of the phase field. For the small frequencies the relaxation part is found to be much 
larger than the usual kinetic part of the phase dynamics which is proportional to the second power of the frequency. 
In order to describe the phase difference relaxation process the consideration is developed in the framework of the 
Schwinger-Keldysh technique [10], [11] which is very convenient for the description of quantum kinetic processes. 

The generating functional of the system can be written in the form Z = ^ D^D^eexp {i {Sjj + Sj)}. The action 
Su of the inhomogeneous non-equilibrium Bose gas is given by 



Su = 



j^dtj dPr jt/- [idt + fi{x)-U {x)] V - i (yi) (^V) - ^ 



and the term Sj reads as Sj = § dt J d^r [i/jj + jif) , where j and j are the infinitely small sources. The integration 

over the time t in the action is realized over the time reversing contour and is denoted as § dt. The Planck constant 
h and mass m of a Bose particle are set equal to unity h = m = 1. The chemical potential jj, (x) is supposed to have 
the constant value /xi = /x — ^A/x for the coordinate x>d and /X2 = /x -f ^A/x for x<-d where the chemical potential 
/X obeys the equality fi = \n and the chemical potential difference obeys the equality A/x = AAn. The small value of 
An results in the small value of the chemical potential difference A/x << /x. 

Due to the large magnitude of potential Uq in the region xG [— d, d] the interaction between Bose particles can be 
neglected and the integral for Z over the fields ipi'^l' takes the Gaussian form in this region of the x-coordinate. In 
this connection the integral for Z over the fields i/; (x, ~r^_L; t) ,ip {x^l^ ^]t) for xG [—d^d] with the fixed magnitudes 
of these fields at the boundary of the region of the nonzero external potential can be calculated. Moreover, the 
characteristic scale of the i/)-field variation in the time and in the 'r' ±_ space is supposed to be much larger than 1/C/o s.nd 
l/VT^, respectively. It is these fields that describe the slow relaxation process possessing the characteristic frequency 
proportional to the small tunneling coefficient. Due to this assumption the tunneling amplitude of the if) fields has 
the local character for the time variable t and space variables ~r' ±. The magnitudes of the i/;-fields at the boundary of 
the region of the nonzero external potential are denoted as ■0 (d, ~r*_L, t) = ^si {'^±jt); ip {—d, l^±,t) = ipa2 {'^±j t). 
The calculation of the integral for Z over the ^'-fields in the region xG [—d, d] gives the action in the form of the sum 
of four terms S = Syoi 

+ Siurf + S^surp^ + ^i- The part of the action 

S,oi = jdt j dxj d^ri_ I {idt + m) ^ - ^ ( (^V') - \ 'I (1) 

C/(x)=0 

where the region of the integration over the x-coordinate is the sum of the regions {—L, —d) U (d, L). The parts of the 
action S'^^^.j and Sf^^p^ can be represented in the form [17] 

S[t} = ~coth{2Hd) j dt j rfV {! V.i I' + I ^s2 (r^,t) |2} (2) 



where the magnitude x is x = \/2Uq, the index 1 denotes the right-hand side of the bulk and the index 2 denotes 
the left-hand side of the bulk. The term Sj can be written as Sj = § dt ^ (V'sai + j'4'sa) where j and j are the 

a=l,2 

infinitely small sources. Later on we assume that the width d of the potential barrier is sufficiency large so that the 
inequality xd >> 1 takes place and, thus, sinh (2>£:d) « cosh(2>f(i) « ^e'^'^'^ » 1. The term S',^"""'' describes the 
tunneling between the right-hand and left-hand sides of the bulk. From the form of this tunneling part of the action 
it can easily be seen that the tunneling amplitude for the condensate particles and for the low energy non-condensate 
particles has the same nonzero magnitude in contrast with the work [8]. 

Our goal is the calculation of the functional integral for Z over the ^/i- fields in the bulk, i.e., in the region xG 
(— L, — d) U (d, L) with the fixed values at the potential barrier boundary, and, thus, obtain the effective action for the 
fields at the boundary ipsi, i's2- 

Below the onedimensional character of the system and its excitations is supposed. In this supposition the considering 
fields depend on the x coordinate only. Later on, it is convenient to introduce the dimensionless space and time 
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coordinates via the following change of these variables x ^x: t ^ t/fi. Thus, the system of units which we use 
later measures the quantities of the dimension of length in units of ^ and the quantities of the dimension of energy 
in units of fj,. Moreover, we translate the x coordinate for the right-hand side of the bulk on the value -d and for the 
left-hand side on +d, so as the right-hand and left-hand boundaries of the potential barrier take the zero coordinates 
x=0. _ 

It is convenient to represent the fields ip in the modulus-phase form ^ = Re"^; ip = Re~^'^. Below the modulus R 
of the field tjj is represented as the sum of the saddle-point configuration F {x) and the fluctuations of the modulus 
field p which characterize the density fluctuations 



R{x,t) = F{x)+p{x,t) 
The saddle-point configuration F (x) obeys the equation 

1_ 
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F^ (x) ]F{x)= kF {x) S (x) 



(4) 



(5) 



and has the form F (x) ~ tanh (| x \ +Xo), and the constant Xq is the small value Xq 
be represented as the sum of two summands 



p {x, t) = Sp {x, t) + 



Vcc-F {x) 
V,F(0) 



Ps{t) 



The fluctuations p can 



(6) 



where the fluctuations 5p {x, t) obey the zero boundary conditions 5p (0, t) — 0. Note, that the second summand of (6) 
is the zero mode of the Hamiltonian Hi = — ^V^ — 1-1- which describe the dynamics of the fluctuations p {x, t). 
To describe the slow process of the tunneling through the potential barrier we are interesting for the flelds 6p and 
which are slow varying values at the space scale ^ and the time scale l/p.- For the derivation of the slow fleld effective 
action we substitute the field ijj into the action (1-3) in the modulus-phase representation where the field R(x,t) has 
the form (4) and where p is determined by (6). At the same time, in certain terms of the action due to the slowness of 
5p and (p the function like F^ {x) can be replaced by the function 1 — 6 (x), where 6 {x) is the 5- function. Integrating 
over the fields 6p {x, t) and ps (t) in the generation functional we obtain the effective action for the slow phase field 



= K jdt 



dx 
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6{x) 



U(x) 



where 7 = ;^ exp {—2>cd), the value Aips (t) is the difference of phases at the right-hand and the left-hand boundaries 
of the potential barrier A(ps (t) = <fs2 (t) — fsi (t), the constant K is K = n^S± » 1 and 5 (x) is the Dirac 5- function. 
The value S± is the square of the bulk section which is supposed to obey the inequality S± < to ensure the ID 
character of the system. The dimensionless constant K = n^S± has the large value K » 1 due to the large value of 
the parameter n^^ >> 1. The last inequality is ensured by the application of the gas approximation. 

The last step to the derivation of the effective action for the phase difference A(fs (t) is the integration in Z over 
the bulk components of the phase (p{x,t) with the fixing boundary values cps For this purpose the field ^p{x,t) 
is represented in the form of the Fourier expansion ip{x,t) = '^<fik {t) e'*^^ with the Fourier components (fik obeying 

fe 

the condition ipg (t) =^<Pk (i) where the coordinate of the right-hand boundary of the potential barrier we choose 
fe 

as x=0. The integration over ipk gives 



Seff [fs 



j> dtdt' 



^-{dt^,f5{t-t') + ]^^s {t) 



<D> 



ifis {t') + 7 cos {A(fs {t) - Aij,t) 6{t- t') 



(8) 



where we denote < D >= J ^Dk- In the frequency-momentum representation in the "triangular" form [10], [11] at 
the temperature T the phonon propagator D {w, k) can be written 



D{Lu,k) 



D^{u,k) 



where D^'^ {u, k) = {u ± iSf - P and (w, k) = coth 
give 



D^{uj,k) D^{u>,k) 

jjR _ jjA\ rpjjg simple calculations of < 5 > 



(9) 
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<D^>^ 



where < 



1 



(10) 



II. THE EFFECTIVE ACTION FOR SMALL FREQUENCIES AND THE FOKKER-PLANCK EQUATION. 

Introducing the half-sum and the difference of the fields instead of the fields ipgi and (pg2 we obtain the effective 
action for the difference of phases at the boundary of the potential barrier. In the "triangle" representation it can be 
written as 



•[p]=2K j 



dtlf 



1. 



-uj + iu! 



3> + $ 



1, 



-ui 



ip + lif 



2a} coth 



2T 



(11) 



where ui = idt, the function V (p] is V <p] = 7 sin ($ + A^.t) sin (<p). In this representation the fields $ and ip 
are expressed via the fields belonging to the upper and the lower time contour branches as $ = i + ip^) and 
<p = ^ {>fi+ — ip-). Note that as usually in the "triangle" representation of the Keldysh-Schwinger technique the field 
$ describes the kinetics of the phase difference and the field (p describes the quantum noise of this value. 

The fields $ and p can be fragmented into fast and slow components. The slow fields are the fields having the 
frequencies u obeying the inequality u> < coq and the fast fields are the fields with the frequencies u> > coq where 
7 << Wo << 1- The part of the action for the fast fields is obtained by the expansion of the potential V up to the 
second order of the fast components 5$, d(p and has the form 



S'^f""'^ [6ip, S^] = 2K 



. Sp [iS^ + iu)] 5^ + (5$ [icD^ - ioj] 6p + iSp [2ui coth (^)] 6p- 1 



(12) 



V^'^ = 7 cos ($ + Ajj^t) cos {p) ; V^^ = V^^ = -7 sin ($ + A/xf) sin {p) 

where Sp, 5^ are the fast components of fields and p, $ are the slow components. The part of the action for the slow 
fields can be written as 



S(siow) =2K J dt l^-2pdt^ - V [$, if] + ip> 



2u) coth ( I 
2T) 



(13) 



Due to the large value of the constant isT >> 1, so as the value ;^ In :^ obeys the inequality ;^ In ;^ << 1, the action 
gifast) approximated as 

[5^, 5^] = 2K [dt{^^ + + [g^ - + *<5^ [2J coth (#)] Sp>- | 

^ ^' ^ J \ -V''^ < 6p6^ >o -^V^^ < >o -iV^;; < S^Sp >o j ^ ^ 



where 



\lj\>uio 



(15) 



\ij}\>u:o 



The correlator Do is defined by the equality Dq = (^^'(Jx + ^ < -D > , where is the a^- Pauli matrix. Thus, 
taking into account the renormalizations resulting from the fast fluctuations the effective action for the slow fields 
can be written as 



Si^iow) ^]=2K dt\ -2p (dt^) - V [$, p] + 



■ ip 



2u) coth I — - 
2T 



(16) 
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In this expression the renormalizations proportional to the correlators < S^tS^t >o give the small renormalization 
of the constant 7 in the potential V. They can be neglected due to the smallness of the value In ^3^^ << 1- On 
the other side the renormalizations proportional to the correlators < Sipt6^t >o are taken into account because they 
can become comparable with the relaxation term ip [2a;coth (^y)] V in the case of the low temperatures and small 
frequencies T ~ , w ~ . 

It is convenient to shift the field $ ^ $ — Aiit. As a result of this shift the action (16) takes the form 



2wcoth I — 



where ^ = 7 sin $ sin [ip) — VQip, at that Vq = 2A/x. Taking into account the smallness of the quantum noise (fi 
due to the large value of the constant K » 1 the action can be expanded in the powers of the field (p. This expansion 
gives 



■ lip 



2Q coth I — 



1 



16K 



7 cos $ 



where 



= 7sin$ - uo; A'($)=7cos$ 



^-Ai<i>)<p-t^A'{^)\ (17) 



(18) 



The first correction to the renormalized amplitude of the potential V (p], i.e., to the constant 7 with respect to the 
value 1/K has the form 7* = 7 ^1 — In ("^^^ • Integrating over (p in the generation functional Z we obtain 



q{sIow) 
^eff 



[$] = -i^pln (f) +2kJ dtdt'i^^[-2 {dt^) - At] [-2 {df^) - At^] - i^A'^S {t - t') 
where the correlator Ft^c is 



r = 



22 coth 



1 



;7* cos$ 



2T J 16K 

In the approximation of the local in time action wc have 



2u) coth 



2T 



1 



16K 



1 



Si^'p^ [$] = i2K JdfJ^\ [2 {dt^t) + A ($,)]' f - ^ A' ($,) I 



(19) 



(20) 



r 1 _ 1 

° 4r+^7* 4T* 

Finally, the action can be represented in the form 



Sifr' m = J dt[^im' + ^A- - ^A' 



(21) 



(22) 



where M = 2To = ■ Note that the generation potential Z with the action •S'^//"' [^] (^2) has the form correspond- 
ing to the quantum mechanics with the "imaginary time" in which the role of the "imaginary time" plays the usual 
real time. This action corresponds to the action which is obtained for the classical random process describing by the 
Langevin equation with the white noise [14] and having the form of the supersymmentric quantum mechanics. The 
Schredinger equation and the Hamiltonian corresponding to the action S'^^j"'^ [$] are 



~2^^** = Hfp^ 



(23) 



Hpp = — 



1 d^ ^a2_J_A' 

{2Kf 2Mrf2$ +8 B,K 



(24) 
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This Hamiltonian Hpp can be represented in the factorizcd form 

Hfp = - 2M{2Kr - MKA) + MKA). The transformation * = exp {Mif^ ($)} P ($), where F ($) = 
7* (1 — cos$) — vq^ and A = (V ($)) transfer the Hamiltonian (18) to the form of the Fokker-Planck operator 
Lfp = — 2M(2K)'^ 'M "^MKA^. The stationary solution of the corresponding Fokker-Planck equation in the case 

of the zero velocity vq = Q has the form P ($) = Co exp ^— ^^i*^ ^ , where V ($) = 7* (1 — cos $) and the constant 

Co is Co = ^ • In the case of the nonzero velocity vq the stationary solution of the Fokker-Planck equation is 

determined by the equation 



(J^+2MKA^P{^)=j 



and can be found taking into account the condition of the periodicity of the probability distribution P (0) = P ($) and 

27r 

the normalization condition J (i$P($) = 1 [13], [18]. The constant j is found from the normalization condition for 



P ($) and for the case of the small An wc obtain j — ^ exp |— 2 (^^^^ |. The probability of the phase slip process is 

determined by the value {dt < $ >) ^ j. For the zero temperature and small velocities vq, so as vq = A/i << T* , the 
value T* can be considered as the effective temperature and for the large time t, when the action can be considered 
in the local in time form, determines the diffusion coefficient. In the opposite case = A/i » T* the diffusion 
coefficient is defined by the value Tt^t = (^^) ^^"^ this value should be substituted instead of the value {"/* /T*). 

This can be seen from the expression (20) for F with taking into account that in the case of the large value of A/x 
compared with T* the frequencies iv should be cut of for the small values by the chemical potential difference A/i. 
Note, that the result for the probability of the phase slip (l/rps) ~ wo exp 1—27^ In | coincides with the result 

of the works [15-17]. 

The Fokker-Planck equation describing the relaxation process of the difference of phases has been obtained both for 
high temperatures and for low temperatures. For high temperatures T" > 7 the obtained Fokker-Planck equation is 
equal to the classical Langevin equation with the white noise with the correlator proportional to 1/T [18]. For the low 
temperatures T << 7 and the small density difference An, so as the inequality vq << 1/K takes place, the relaxation 
kinetics of the difference of phases is described by the Fokker-Planck equation with the effective temperature T* and 
the diffusion coefficient determined by both the usual temperature and the quantum fluctuations of the Bose gas. For 
the small temperatures and the large value of the density difference An, so as iio >> 7/i^, the diffusion coefficient is 
determined by the value vo- In the last case the result is analogues to the result of the works [15-17]. 

The author thanks S. Burmistrov, A. Kozlov and Yu. Kagan for helphil discussions. This work was supported by 
the Russian Foundation for Basic Research, by the Netherlands Organization for Scientific Research (NWO) and by 
INTAS -2001-2344. 
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